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• Abstract 

CN ' Let V = A^V be the iV-fermion Hilbert space with Af -dimensional single particle 

' space V and 2N < M. We refer to the unitary group G of as the local unitary 

• . (LU) group. We fix an orthonormal (o.n.) basis \vi), . . . , \ vm) of V . Then the Slater 

Q I determinants ej^^..._ijy := \vi^ Avi^ A • • • At)j^) with ii < • • • < form an o.n. basis of V. 

I Let 5 C V be the subspace spanned by all ei-^^,,,^if^ such that the set {ii, . . . , ijy} contains 

^ ■ no pair {2k — 1, 2k}, k an integer. We say that the IV') G S are single occupancy states 

, (with respect to the basis \vi), . . . , \ vm))- We prove that for = 3 the subspace S is 

^ I universal, i.e., each G-orbit in V meets 5, and that this is false for > 3. If M is even, 

^ ■ the well known BCS states are not LU-equivalent to any single occupancy state. Our 

- - - main result is that for = 3 and M even there is a universal subspace W C 5 spanned 

by M{M - 1)(M - 5)/6 states ei^^,„^ij,. Moreover the number M{M - 1)(M - 5)/6 is 
minimal. 



1 Introduction 

Entanglement is at the heart of quantum information theory. It makes possible secure and 
high rate information transmission, fast computational solution of certain important prob- 
lems, and efficient physical simulation of quantum phenomena [H [2] . The most fundamental 
property for any kind of study of entanglement is that it is invariant under local unitary (LU) 
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transformations. The celebrated Schmidt decomposition for bipartite pure states provides a 
canonical form for these states under LU, which makes possible a complete understanding 
of their entanglement properties. In particular, the entanglement measure is characterized 
by the entropy of Schmidt coefficients. 

For multipartite systems, however, no direct generalization is possible. The study of 
entanglement for multipartite pure states is hence much more challenging than the bipartite 
case. Considerable efforts have been taken during the past decade. However no complete 
satisfactory theory can be reached as for the bipartite case [2|,[3]. Still, the first step toward 
the understanding of multipartite entanglement is to study the orbits of pure quantum states 
under LU. Acin et al started from looking at the simplest nontrivial case of a three qubit 
system [U |5]. They introduced the concept of local basis product states (LBPS), which are 
three-qubit product states from a fixed set of single particle orthonormal basis. They studied 
the property of subspaces spanned by LBPS such that every three-qubit pure state is LU 
equivalent to some states in such a subspace. This kind of subspace can hence be called 
'universal subspace' [HE]. 

Acin et al. showed that five LBPS are enough to span a universal subspace for three 
qubits [1]. For one such choice given by {|000), |100), |101), |110), |111)}, they can further 
restrict the regions for the superposition coefficients such that it gives a canonical form 
for three qubits [6]. That is, each LU orbit corresponds to only one set of parameters in 
the canonical form. Generalizations to > 3 qubits are developed and it is found that 
a universal subspace needs to be spanned by 2^ — A^ LBPS, and such a minimal set of 
LBPS are identified [7]. Unfortunately, this does not provide much saving compared to the 
2^-dimensional Hilbert space. This is simply due to the fact the LU group is such a small 
group compared to the entire unitary group U(2^). 

Entanglement theory for identical particle systems has also been considerably developed 
during the past decade [3], [8], |9l [TOl [TTl [12]. Entanglement for identical particle systems 
cannot be discussed in usual way as for the distinguishable particle case. This is because 
the symmetrization for bosonic systems and antisymmetrization for fermionic systems for 
the wave functions may introduce 'pseudo entanglement' which is not accessible in practice. 
It is now widely agreed that non-entangled states are reasonably corresponding to the form 
|t>)^ for bosonic states [HI [3l [12] (there was indeed suggestion that non-entangled bosonic 
state corresponds to \vi) V \v2) V ■ ■ ■ V \vn) [10], however most of the literatures accept \v)'^) 
and to the form of Slater determinants \vi) A \v2) A ■ ■ ■ A \vn) for fermionic states [31 [12]. 

For bipartite fermionic pure states, a direct generalization of Schmidt decomposition 
is available [9]. That is, every bipartite fermionic pure state is LU equivalent to a Slater 
decomposition ^ lA), where = (Al/^i) = ^ij and {ai\(3j) = 0. This thus 

allows a complete understanding of entanglement properties for bipartite fermionic pure 
states. Many interesting analogues to the bipartite distinguishable particle case have also 
been identified, such as concurrence and magic basis [HI [12] . The entanglement of formation 
(EOF) for mixed states has also been investigated for the single particle state Hilbert space 
of dimension four. In particular, a formula was identified [9l [T2] and is similar to the known 
formula for two-qubit states by Wootters [T3] . 
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The generalization to more than two fermions has also been discussed [12]. However 
similar to the distinguishable particle case, there is no Slater decomposition for multipartite 
fermioinic system. Then it is more difficult to study their entanglement properties. The ffist 
step is still using the LU orbits. Due to particle-hole duality, in most cases one will need 
to have M > 2N to obtain meaningful discussions, where M is the dimension for the single 
particle state Hilbert space V. For examples with M > 2N, see the BCS states defined in Eq. 
flH|) . The reason is that the Pauli exclusion principle requires M > N, and the particle- hole 
duality gives equivalence between N, M and M — N, M systems. Also, it is even nontrivial 
to see whether a given fermionic pure state is unentangled (i.e. LU equivalent to a single 
Slater determinant), where one needs to check the Grassmann-Pliicker relations (see e.g. p. 
A 111.172 Eq. (84-(J,H)) in [H], Prop 11-32 in [15], and [I2]). 

In this paper, we take one further step beyond just the Grassmann- Pliicker relations for 
the entanglement properties of multipartite fermionic pure states. We study the universal 
subspaces spanned by Slater determinants built from orthonormal local basis. So we can 
build a fermionic analogue of the universal subspace spanned by LBPS in the distinguishable 
particle case. The obvious difficulty is the above-mentioned M > 2N condition, for which 
one can no longer hope for any elementary method (e.g. linear algebra) as was useful in the 
iV-qubit case. However, this difficulty, once overcome, might become an advantage given 
the relatively large LU group U(M) compared to the entire unitary group on the A^-particle 
fermion space V = A^V which is of dimension (^). In other words, one can ideally hope 
for a 'saving' of the dimension of U(M), which is of the order M^. That is, the best one can 
hope for the dimension of a universal subspace is of the order (^) — for N > 2. 

Our focus will be mainly on the N = 3 case, where we show that such a saving of order 
is indeed achievable. Not only we want to find a universal subspace which mathemat- 
ically achieves the order (^) — M^, but also we want the universal subspace with a clear 
physical meaning. To obtain such a universal subspace, we introduce a configuration for the 
Slater determinants called single occupancy. This is typical in condensed matter physics for 
studying properties of strongly-correlated electron systems, for instance high temperature 
superconducting ground state of doped Mott insulator (see e.g. [TB]). 

Single occupancy fermionic states are also of interest in quantum information theory 
in recent years, as there exists a qubit to fermion mapping which is one to one between 
qubit states and single occupancy fermionic states. This mapping was used to show the 
QMA completeness of the fermionic N-representability problem, that is, the overlapping N- 
representability problem is hard even with the existence of a quantum computer (TTJ [50] . 
The same mapping was also used to disprove a 40-year old conjecture in quantum chemistry 
using methods of quantum error-correcting codes [T8] . 

The physics considerations lead to the following picture. We illustrate this in Fig. [TJ 
Suppose we have in total K sites. At each site we have three possible states: 1) no spin (a 
hole, shown as a dot); 2) spin up (shown as arrow up); 3) spin down (shown as arrow down). 
Now suppose we have in total spins. Then there are in total (^^) spin configurations 
(note that once the spin states are fixed, the holes are fixed). As an example. Fig. [2] shows 
the case A^ = 3 and K = 4, i.e. 3 spins in 8-dimensional space. Figure A is a single 
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occupancy configuration which has at most 1 spin on each site. A single occupancy state 
is a superposition of single occupancy configurations. Figure B is not a single occupancy 
configuration where site 1 is occupied by 2 spins (i.e. double occupancy on site 1). While 
the single occupancy configuration is defined for M = 2K which is even, one can readily 
generalize the configuration to the case of odd M, leaving one single-particle spin state 
unpaired. 




1 2 3 K-\ K 

Figure 1: Physics picture for single and double occupancy 




1 2 3 4 1 2 3 4 

A B 

Figure 2: 4 sites 



To be more precise, for the single particle Hilbert space V of dimension M = 2K, which is 
a tensor product of the spatial part with dimension K (with a fixed basis \gi), \g2), . . . , Igx), 
which each \gj) highly localized on the jth site) and the spin part with dimension 2 (with a 
fixed basis | t); I i)), then a single occupancy state is a Slater determinant where \gj) ® | t 
)) \9j) ® I i) do not show up at the same time for any j = 1,2, . . . K. For instance, the single 
occupancy state shown in Fig.[2JA. can then be written as (1(71)®! i))/\i\g2)®\ t))/\(lfi'4)® I i))- 
However, the state given in Fig. [2f3 is {\gi) ®\t))^ ildi) ® | i)) A {\g2) ® | t))) which is not 
single occupancy as both l^fi) ® | t); Ifi'i) ® | i) show up. 

For a simplified notation, we can relabel the 2K single particle basis {Is'j)®] t); Ifi'i)® I i)} 
as {\2j — 1), |2j)}, for j = 1,2 . . . , K . The M odd case is hence similar, just that one can 
correspond the first M — 1 basis states, denoted by {|2j — 1), |2j)} to {\gj) ® | t); \9j) ® I i)} 
for j = 1, 2 . . . , (M — l)/2, and then the last 'unpaired' basis state |M) to |5'(m+i)/2) ® | t); 
as the state \g{M+i)/2) ® | i) is 'inaccessible' by any of the A^-fermions. Note that despite 
the choice of the labelling, the basis of the Hilbert space is fixed throughout the paper. This 
corresponds to a simple fact that the spatial sites are physically fixed. 

Our main theorem of the paper proves that for = 3, the single occupancy subspace 
S (i.e. subspace of all single occupancy states) is universal. In other words, any A^ = 3 
fermionic state is LU equivalent to a single occupancy state. Here by LU we mean that local 
unitary transformation on the entire M-dimensional single particle Hilbert, i.e. elements in 
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U(M), not only unitaries on the spatial part (or the spin part) of the wavefunction only. This 
definition of LU agrees with those used in quantum information community for fermionic 
systems (see, e.g. [SI [12]), which preserves entanglement properties, such as the spectra of 
reduced density matrices. 

For M even, we can also obtain universal subspaces contained in S whose dimension 
is minimal. Our main tool is to apply Theorem 4.2 of [6] to the case of LU groups with 
representation on the iV-fermion Hilbert space V = A^V^. It should be emphasized that 
this application in our case is highly non-trivial. Indeed, after selecting a candidate for a 
universal subspace, one has to prove that the main condition of the theorem is satisfied. For 
small M this can be done by using a computer, but the proof that it can be applied for all 
M is hard to find. Finding a universal subspace of minimum dimension is much harder and 
we succeeded to find one which works for even M. For the case when M is odd, one may 
conjecture that the candidate subspace described in Proposition [S] (see Eq. is universal 

but this has been verified only for M = 7, 9, 11. 

For the smallest non-trivial case of = 3, M = 6 where DimA'^I^ = 20, a direct 
application of our results gives universal subspaces with dimension 5. This is indeed a 
big saving. It is interesting to compare this result to Acin et. al's result of three qubits, 
where they also have universal subspaces of dimension 5, which are spanned by LBPS. This 
may be related to the fact that, when M = 2N, each qubit state corresponds to a single 
occupancy state based on the qubit to fermion mapping mentioned above. It should also 
be mentioned that, although it is not directly applicable, our work may shed light on the 
understanding of the N-representability problem [19] , where the single particle eigenvalues 
are invariant under LU. Indeed, our universal subspace with dimension 5 for the case = 
3, M = 6 gives an alternative proof of the N-representability equalities and inequalities in 
that case [20l|2ll [22l|23] . 

To complete our study of single occupancy, we further prove that the single occupancy 
subspace is not universal for A^ > 3, i.e. not all fermionic states are LU equivalent to a single 
occupancy state, though for fixed A^ almost all states are single occupancy in the large M 
limit. Our argument is based on dimension counting. For concrete examples, we show that 
for A^ even, BCS states are not LU equivalent to any single occupancy state. This is intuitive 
as BCS states are always 'paired' so they are most unlikely to be transformed with LU to 
something unpaired. 

We organize our paper as follows. In Sec. [21 we provide preliminaries used throughout the 
paper. We formally define antisymmetric tensors, decomposable states (Slater determinants), 
LU operations, universal subspaces and single occupancy states. We also discuss the lower 
bound for dimension of universal subspaces. It sets the goal for later work. In Sec. [3l we 
consider universal subspaces for A^ = 3. We prove our main theorem that for A^ = 3, the 
single occupancy subspace S is universal, see Theorem [H For M even, we can also obtain 
subspaces of S whose dimension is minimal, and explicit choice of Slater determinants for 
those subspaces are given. In Sec. [H we try to generalize our theorem to the A^ > 3 case 
but obtain a negative result. That is, we prove that the single occupancy subspace is not 
universal for any M > 2N. So there always exist some states which are not LU equivalent 
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to single occupancy states. For even, we use BCS states as concrete examples of states 
that are not LU equivalent to any single occupancy states. Finally, a brief summary and 
discussion will be given in Sec. O 

2 Preliminaries 

2.1 Decomposable states and inner product 

Let y be a complex Hilbert space of dimension M and % = ®^V^ the A^th tensor power of 
V. The inner product on V extends to one on T-L such that 

N 

{Vi, V2,..., Vn\Wi, W2,..., Wn) = Yl{Vi\Wi). 

i=l 

We denote by A^\^ the A^th exterior power of V , i.e., the subspace of "H consisting of the 
antisymmetric tensors. We refer to vectors \ip) G /\^V as N-vectors. We shall often identify 
the A^- vectors with antisymmetric tensors by using the embedding A^V — ?• "H given by: 

\vi) A ■ ■ ■ A\vn) ^ ^ sgn(cr)|t;^(i),'y^(2), . . . ,u^(Ar)), (1) 
o-eS'jv 

where is the symmetric group on A^ letters and sgn[a) the sign of the permutation a. 
We say that an A^-vector {ip) is decomposable if it can be written as 

= \vi) A--- A |^;7v) 

for some \vi) G V. For brevity we may also write 

I^A) = \viA---Avn) 

when no confusion arises. Decomposable states are also called Slater determinant states or 
Slater determinants in physics. For the sake of brevity we shall mainly use the name of 
decomposable states. Unless stated otherwise, the states will not be normalized for conve- 
nience. 

If is a vector subspace of V, then A^W is a vector subspace of A'^V. Given a 
\ip) G A^V, there exists the smallest subspace W <^V such that \ip) G A'^W. We shall refer 
to this subspace as the support of \ip). By Eq. ([T]), this is equal to the range of the reduced 
density matrix of any single particle of {tp). In the case when {tp) 7^ is decomposable, 
say \ip) = A ■ ■ ■ A V]y), then its support is the subspace spanned by the vectors \vi), 
1 = 1,. ..,N. 

In general, an A^-vector can always be written (not uniquely) as a sum of decomposable 
A^-vectors. With this notation, the inner product in a'^V is characterized by the fact that the 
inner product of two decomposable A^-vectors \ip) = \wi A ■ ■ ■ Awn) and \ip) = \vi A ■ ■ ■ A Vn) 
is equal to the determinant of the N x N matrix [(mjltij)]. It is also equal to -^{ipltp) where 
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the inner product is computed in the tensor space by using the identification given by Eq. 
([T]). Hence, \ip) _L when this determinant is zero. 

The 2-vectors G /\^V are often identified with antisymmetric matrices of order M 
via the isomorphism which assigns to any decomposable 2- vector = \v A w) the ma- 
trix \v){\w)'^) — \w){\v)'^). Under this isomorphism, the nonzero decomposable 2-vectors 
correspond to antisymmetric matrices of rank 2. 

We shall use the well known exterior multiplication map from A"'V x A^V into A'^^"'V, 
which is a complex bilinear map. Note that this map is identically zero ii N + n > M. Let us 
also introduce the 'partial inner product' operations. Given an n-vector \ip) and an A^-vector 
lip) with n < N, we can define their partial inner product {(f\ip) G A^~"V. In other words, 
this is a map 

A^'V X A^l^ ^ A^-^'V 

which is antilinear in the first argument \(p) and complex linear in the second argument 
Thus, it suffices to write the definition in the special case when \if) and are decomposable, 
say \ip) = \wi A ■ ■ ■ A w„) and = \vi A ■ ■ ■ A v^). In that case we have 

= ^sgn{a) (^{Wi\v^i)^ \v^^n+i) ^- ■ ■ Av^^N)), (2) 

where the summation is over all permutations a ^ Sjy such that the sequence a(n + 
1), . . . ,cr[N) is increasing. This definition is a modification of the one given in [TH p. A 
III. 166] for the partial bilinear inner product. We have omitted the overall factor (— 
and replaced the bilinear inner product with the Hilbert space inner product. More impor- 
tantly, it appears that there is an error (or misprint) in Bourbaki because the right hand side 
of their formula is not antisymmetric in the arguments x*. We had to modify their formula 
by requiring that only the second one of their two sequences of a-values is increasing. The 
Eq. ([2]) can be rewritten as 

{ip\lp) = '^SgIl{a){wiA---AWn\v^(i)A---Av^^ri))\Va{n+l)/\---/\Va{N)), (3) 
a 

where the summation now is over the permutations cr e Sn such that both sequences 
cr(l), . . . , a{n) and cr(n + 1), . . . , cr{N) are increasing. 
If n = 1 then the formula reads as follows: 

N 

{wi\ip) = ^{-iy~^{wi\v,)\viA---AVi.iAVi+iA---AvN). (4) 

i=l 

If we define the interior product l{w) to be the linear map A^^^V which sends 

to {w\ip), then one can check that Eq. ([2]) can be rewritten as follows: 

(V?!^) = L{wn) o ■ ■ ■ o l{wi){\vi A ■ ■ ■ Avn)). (5) 

By using the embedding ([T]), the partial inner product given by Eq. ([2]) is a positive scalar 
multiple of the inner product {(plip) where both \ip) and are viewed as antisymmetric 
tensors. Note that in the case n = N, the partial inner product in Eq. ([2]) is the same as 
the inner product on a'^V defined at the beginning of this section. 
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2.2 Local unitary equivalence 

U Ai : V ^ V , i = 1, . . . , N , are linear operators, then their tensor product ^fLiAi will be 
identified with the unique linear operator on "H which maps 

\vi) ® \v2) ® ■ • ■ ® \vn) Ai\vi) ® A2\v2) ® ■ ■ ■ ® An\vn), \vi) eV, (z = 1, 2, . . . , N). 

li Ai = A2 = ■ ■ ■ = An = A, then we shall write ^^A or instead of ®^iv4j and refer 
to as the iVth tensor power of A. U \tp) G A^V then we also have ^'^A{\tp)) e A'^V. 

Consequently, we can restrict the operator ®^A to obtain a linear operator on A^^, which 
we denote by A'^A or A^'^ . We refer to A'^ A as the iVth exterior power of A. Explicitly, we 
have 

A^A{\vi A t;2 A ■ ■ ■ A Vn)) = \Avi A Av2 A ■ ■ ■ A Avn)- 
The general linear group G := GL{V) acts on by the so called diagonal action: 

A-{\vi)^\v2)^---^\vn)) = A\vi)cqA\v2)®---^A\vn), AeG, \vi) e V, {i = 1,2, . . . , N). 

In other words, A & G acts on "H as ®^A. Similarly, we have the action of G on A^^ where 
A & G acts as A'^A. For convenience we shall abbreviate both of these two actions by a 
i.e., we have 

A ■ \vi, . . . ,vn) = \Avi, . . . ,Avn), (6) 
A- \vi A ■ ■ ■ Avn) = \AviA--- AAvn)- (7) 

Both of these actions can be restricted to the unitary group U(y) of V. 

We shall say that two A^- vectors |0) and {tp) are equivalent if they belong to the same G- 
orbit, i.e., = A - |0) for some A ^ G. We shall also say that they are unitarily equivalent 
or LU-equivalent if such A can be chosen to be unitary. For example, any decomposable 
iV-vector is unitarily equivalent to a scalar multiple of |1) A ■ ■ ■ A |A^). 

In practice, LU can be realized by a Hamiltonian 

N 

H = Y,H„ (8) 

i=i 

where Hj acts on the jth particle, Hj\vi, . . . , vjy) = \vi, . . . , HjVj, . . . , v^)- The evolution of 
the system is unitary, i.e. 

N N 

e-^^* = exp {J2 -^H,t) = (g) f/„ (9) 
j=i i=i 

with Uj = e"*^-**. The Hamiltonian of the form Eq. ([H]) is usually called a single-particle 
Hamiltonian. For a fermionic system, we have Hi = H2 = ■ ■ ■ = = H , hence Ui = U2 = 
■ ■ ■ = Un = U. Hence the evolution of the system is indeed the fermionic LU A^U. 
In the case = 2, we have the canonical form for unitary equivalence. 
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Lemma 1 If N = 2, then any 2-vector \ip) is unitarily equivalent to Yli=i '^d^^ — 1) A \2i) for 
some Ci > ■ ■ ■ > Cfc > and some integer k > 0. Moreover, the coefficients Ci are uniquely 
determined by 

Proof. The first assertion follows easily from the antisymmetric version of Takagi's the- 
orem, see p4l p. 217]. To prove the second assertion, we denote by K the antisymmetric 
matrix corresponding to {tp). If [/ is a unitary operator on V, then the 2-vector Ulijj) is 
represented by the antisymmetric matrix UKU^ . It follows that the Hermitian matrix KK* 
is transformed to UKK*U\ and so its eigenvalues are independent of U. Since these eigen- 
values are — c^, — , . . . , — c^, — c^, 0, . . . , 0, the second assertion follows. □ 

It can be seen that the Slater decomposition for 2-vectors is essentially the same as the 
well-known Schmidt decomposition for bipartite pure states. While the later can be easily 
generalized to multipartite systems by regarding them as bipartite systems, the same method 
does not work for the former, because the two parties of the bipartite system would have 
different dimensions and this system cannot be antisymmetric. 

We remark that this lemma provides a complete classification of U(y)-orbits on A^V. Let 
us also mention that the antisymmetric version of Takagi's theorem used in the above proof 
was reproved recently in [HI Lemma 1]. However the authors missed to observe that this 
result can be used to construct the canonical form for bipartite fermionic states. Similarly, 
the symmetric version of Takagi's theorem provides a canonical form for pure states of 
two bosons, which is rediscovered recently in [TOl ITTj . See also [25] where more advanced 
mathematical tools were used to solve this classification problem as well as its analogue for 
two bosons. 

In the case when V is 5-dimensional, there is also a simple canonical form for 3-vectors 
under unitary equivalence. 

Corollary 2 If DimV = 5 then every 3-vector {ip) is unitarily equivalent to (ci|l) A |2) -|- 
C2I3) A 1 4)) A 1 5) for some Ci > C2 > 0. Moreover, the coefficients ci and 02 are uniquely 
determined by lip). 

Proof. It is well known that we can write {ip) = \ip) A \v) for some 2-vector \ip) and some 
\v) G V. Clearly we may assume that \\v\\ = 1 and that the support of \(f) is orthogonal to \v). 
Hence, by applying LemmalU we may assume that \v) = |5) and \ip) = ci|l) A |2) +C2|3) A |4), 
Ci > C2 > 0. It remains to prove the uniqueness assertion. Note that \ip) is decomposable 
if and only if C2 = 0, in which case we have Ci = \\ip\\. So, we may assume that C2 > 0. 
Assume that we also have U -lip) = I'lp') = \ A |5), where \ip') = c[\l) A |2) +C2I3) A |4) and 
c'l > c'2 > 0. Since |5) spans the kernel of the linear map V — A'^V sending |x) — )■ \tp) A \x) 
and the same is true for \ip'), we deduce that U\5) = z\5) where z is some phase factor 
(see [151 Proposition 11.28]). We deduce that \ip) and \ip') are unitarily equivalent, and the 
uniqueness assertion follows from Lemma [T] □ 
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2.3 Universal spaces and single occupancy states 

Let G be a connected compact Lie group and V a real G-module. In other words, V is a 
finite-dimensional real vector space and we have a linear representation p of G on V. For 
a vector subspace W C we say that W is universal if every G-orbit in V meets W in 
at least one point. As an example consider the conjugation action of the unitary group 
U{d) on the space oi d x d complex matrices: {A,X) — > AXA~^ where A G l]{d). Then 
the well known Schur's triangularization theorem can be simply expressed by saying that 
the subspace of upper triangular matrices is universal for this action. There are many other 
interesting examples of universal subspaces that occur in mathematics and physics. The first 
question one may ask about the universal subspaces is: what is the minimum dimension of 
a universal subspace? We can give a simple lower bound for this dimension. Denote by Gyy 
the stabilizer of W in G, i.e., the subgroup of G consisting of all elements g & G such that 
p{g) maps W onto itself. Then we have the following lower bound (see O Lemma 4.1]). If 
W is a universal subspace then 

DimW > DimV-DiniG/Gw, (10) 

where the dimension of the quotient space G/Gy^ is equal to DimG — DimGw. (All the 
dimensions here are to be interpreted as the dimensions of real manifolds.) A case of special 
interest is when G^; contains a maximal torus, say T, of G. Then G/T has even dimension, 
say 2m, and we obtain the inequality Dim W > Dim V — 2m. Assume now that in fact V is 
a complex vector space of dimension d, and so of real dimension 2c?, while W is still assumed 
to be just a real vector subspace of V. If W is universal, we must have Dim W > 2{d — m). 

In our application we shall take V = A^V^, where V is a. complex Hilbert space of 
dimension M, and so d = (^). We also specify that G = U(M) and so we have m = (^). 
In particular, for = 3, the above inequality becomes Dim^ W > M{M — 1)(M — 5)/3. If 
W is also a complex subspace, we obtain that Dime ^ M{M — 1)(M — 5)/6. 

In order to define the important notion of single occupancy states in a fermionic system, 
we partition the set of positive integers into pairs {2i — 1, 2i}, i = 1,2, . . .. We shall refer to 
these pairs as standard pairs. For any positive integer i we denote by i the standard pair to 
which i belongs. Thus, we have i = + when i is odd and i = {i — l,i} when i is even. 
We shall say that a finite sequence of positive integers (21,^2, . . . ,'ifc) is a single occupancy 
sequence if the standard pairs 11,12, ■■■ ,ik are distinct. 

Next we fix an o.n. basis {\i) : i = 1, . . . , M} of the Hilbert space V . As explained below 
Fig [21 each of our standard pairs i associated with basis vectors \i + l)} for i odd, then 
correspond to {\gj) ® \ t), \gj) ® | i)}, where j = {i + l)/2. Similar for i even, the standard 
pair i associated with basis vectors {\i — 1), \i)}, then correspond to {\gj) ® | t)) \gj) ® I i)}, 
where j = i/2. 

The ('5^) decomposable vectors \ii A - ■ ■ Ain), 1 < H < •■■ < iN < M, form an o.n. basis 
of A^V. Hence, the dimension of the space A^V is (^). If the sequence {ii,i2, ■ ■ ■ ,iN) is a 
single occupancy sequence then we shall say that |ii A ■ ■ ■ A iAr) is a basic single occupancy 
N -vector (BSOV). Note that if M is odd then M = {M, M + 1} and so there is no state in V 
that is paired with \M). Finally, we say that an A^- vector {ip) G A^^ is a single occupancy 
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vector (SOV), relative to the above o.n. basis, if it is a linear combination of the BSOVs. 
Thus, the set of all SOVs is a vector subspace (the SOV subspace) of and the set of 
BSOVs is an o.n. basis of this subspace. Clearly, this definition depends on the choice of the 
o.n. basis of V . A simple counting shows that the number of BSOVs is 2^(j^) if M = 2K is 
even and 2^(^) + 2^-\/_^) ii M = 2K + 1 is odd. For the case N = K and M = 2K, if 
one ensures single occupancy, then there are indeed only two possible states per site (spin up 
and spin down). In this case, any configuration has a natural correspondence to an A^-qubit 
state. 

Lemma 3 For N = 2, M > 3, any state is unitarily equivalent to an SOV. 

Proof. This is a direct corollary of Lemma [H □ 

3 Universal subspaces for = 3 

3.1 The SOV-subspace is universal for = 3 

The main result of this section is that any tripartite antisymmetric state of dimension M > 5 
is unitarily equivalent to a single occupancy state. Recall we say that a vector subspace 
W C A^V is universal if every {ip) G A^V is unitarily equivalent to some \if) G W. We 
denote the unitary group of V by V{V). If an o.n. basis of V is fixed, we shall identify ViV) 
with the group U(M) of unitary matrices of order M. 

Theorem 4 Let V be a complex Hilbert space of dimension M > 5, {|z) : z = 1, 2, . . . , M} 
an o.n. basis of V , and let Cijk = \i) A \j) A \k) . Then the SOV-subspace is universal, i.e., 
any 3-vector G A'^V^ is unitarily equivalent to an SOV. 

Proof. The case M = 5 follows from Corollary [2l So, we shall assume that M > 6. We 
shall write M = 2K when M is even and M = 2K + 1 when M is odd. 

We apply [6l Theorem 4.2] to the problem at hand. Let us first explain what this 
theorem asserts in this concrete case. We consider the representation p of U(M) on the space 
V := A^V . We denote by T the maximal torus of U(M) consisting of the diagonal matrices. 
Thus, if X G T then x = diag(^i, . . . , C,m) where each C,i is a complex number of unit modulus. 
Each basis vector e^k is an eigenvector of T. Indeed, we have x-Cijk := p{x){eijk) = ^iCj^k^ijk- 

A character, x, of T is a continuous homomorphism into the circle group 5"^ := {z G 
C : = 1}. They form an abelian group under multiplication, but it is convenient to 
use the additive notation. This means that if x' ^-nd x" ^''^^ characters of T, then their sum 
X = x' + x" is defined by x{^) = x'{^)x"{^)- The map Xi '■ T ^ S^, defined by Xi{^) = 
is obviously a character of T. The group of all characters of T is the free abelian group with 
basis {xi : 1 < < M}. By using this notation, we have x ■ Cijk = {xi + Xj + Xk){x)eijk 
and we say that the character Xijk '■= Xi + Xj + Xk is the weight of the eigenvector Cijk- An 
eigenvector which is fixed by T would have weight 0, but in V there are no such eigenvectors. 

We embed the character group of T into the polynomial ring V := R[xi, . . . ,xm] by 
sending Xi Xi for each i. Let X be the ideal of R[xi, . . . , xm] generated by the elementary 
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symmetric functions di, . . . , of the Xj. Let B be the o.n. basis of V consisting of all eijk. 
Denote by Bg the subset of B consisting of all BSOVs. Let W be a subspace of V spanned by 
a subset B' C B. We define its characteristic polynomial fu to be the product of all linear 
polynomials Xi + Xj +Xk taken over all triples k) such that Cijk ^ B'. Then [6l Theorem 
4.2] asserts that the subspace U is universal if fu ^ L 

To prove (a) we shall take B' = Bg. Consequently, we have U equals the SOV-subspace 
and the polynomial Qo '■= fu is given explicitly by the formula 

K 

10 = Y\. Y\. {X2i-1 + X2i+ Xj). (11) 
1=1 jj^2i-l,2i 

We shall need another important fact concerning the ideal X. First observe that X is 
a homogeneous ideal, i.e., it is the sum of the intersections Xd := X n Vd, where Vd is the 
subspace of V consisting of all homogeneous polynomials of degree d. We introduce an inner 
product (-I-) on V by declaring that the basis of V consisting of all monomials is an o.n. 
basis. In particular, -L if 7^ e. We denote by Sm the symmetric group on M letters 
which permutes the M variables Xj and point out that it preserves the above inner product 
on V. The following polynomial of degree 6 := M{M — 1)/2 (with the well known expansion) 

p:= Yl (^i - ^i) = sgn(a)x°ia;^2a;a3 " " " ^^7^ (12) 

will play an important role. Namely, for any / G we know that / G X if and only if 
(/b) = 0; see [6l Theorem 4.2]. For f,g&V, we shall write f = g if g — f&X and in that 
case we say that / and g are congruent modulo X. Thus, ii f = g then = {glv)- Note 
that the degree of go is strictly less than 6. For that reason we shall introduce the polynomial 
q = fiqo of degree 6, where /i = XiX-s ■ ■ ■ xm-3Xm-i if M is even and /i = {xix^ ■ ■ ■XM-2Y if 
M is odd. To prove (b) we shall prove that q ^X (which implies that go ^ ^T). 
From the formula fill I) , we obtain that 

K M-2 

% = nE(^2.-l+X2.)''^'"Vf, (13) 
i=l j=0 

where aj'^ is the jth elementary symmetric function of the variables x^ with k 2i — l,2i. 
(By convention, ctq*^ = 1 and a^^\ = 0.) By using the obvious recurrence formula aj = 
crf^ + {x2i-i + X2i)a'^^li + a;2i-iX2jCr]*i2 = 0, we obtain that 

= J = 0,1,..., M-2. (14) 

It follows that 

K /M-2 j \ 

i=l \j=0 k=0 J 
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We shall need the expansion 

A/-2 j M-2 

^(-iy(x + y)^-^-^X^'"V = (16) 

j=0 k=0 k=0 

where x and y are commuting independent variables. Note that a^!^^ = a^M~2~k for each k, 

and ai""^ = c^l, = 1 for M even while aj^^^ = a^^fl, = for M odd. For convenience, we 

set a^^i = 0. We now distinguish two cases. 

Case 1: M is even. By multiplying Eq. (fT3|) by /i and by using Eqs. ( !T4|) and (fT6|) . we 
obtain that 

K M-2 

^ - n E «f ^-^-'4- (17) 

j=l fc=0 

By setting Si^k = ci^k'^ + o!'l^ili-k^2i-i^2i~^~^ fo^ ^ = 0, 1, . . . , — 1, we obtain 

that 

K K 

q = Wisifl + H h Si,ii-_i) = E n (^^) 

i=i / j=i 

where / runs through all functions mapping the set {0, 1, . . . , — 1} into itself. Let us denote 
by Qf the summand corresponding to /. If / is not a permutation, then each monomial that 
occurs in the expansion of qj will contain two variables with the same exponent, and so 
the inner product vanishes. If two of such functions, say / and g are permutations 

then there is an even permutation of the variables which sends qj to qg. We deduce that 
{qf\p) = {qg\p)^ so {q\p) = K\{qid\p) where id is the identity permutation. 

If the variables X2i-i and X2i do not occur in a polynomial r G P, then by using the trans- 
position which interchanges X2i-i and X2i, we obtain that (a;2f-i^~^^2j^b) — ~(^2i-i^2f~^~^'^b); 
and so 

{s^,kr\p) = {a^\^k - «f ^)(4-i^^-^-'rb). (19) 
Since gid = S10S21S32 ■ ■ ■ sk,k-i, we obtain that 

{q\p) = Klli{a[^l,^,-an- (20) 

i=0 

By Lemma [5] below, we have {q\p) 7^ 0. 

Case 2: M is odd. The proof follows closely the one for Case 1 and we shall only briefly 
sketch the main steps. Recall that in this = 0. By multiplying Eq. (|T3l) 

ve obtain th; 

M-2 \ K M-2 

W ^M-2-kk \ _ TT \^ (M) M-kk 



by the new ^ and by using Eqs. f lT^ and (|T6l) . we obtain that 

K M-2 

n E «r 4-t4- (21) 




^k ^2i-l ^2i I 
k=l ) i=\ k=l 
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By setting = ^x^.^^ + a22,x^,_ia:^-'= for k = 1,2, . . . , K , we have 

K K 

^ = n("U + <2 +■■■+ s[,k) = e n (22) 
i=i f i=i 

where / runs through all functions mapping the set {1, 2, ... , K} into itself. Let us denote 
by qf the summand corresponding to /. If / is not a permutation, then {qf\p) = 0. When 
/ is a permutation, this inner product is independent of /. Thus, {q\p) = K\{qid\p) where 
Qid = ■s'ii'S22'533 ' ' ' s'k,k- usiug a similar argument as in Case 1, we obtain that 

m = Klfliairi.-ar)- (23) 

1=1 

By Lemma [5], {q\p) 7^ which completes the proof. □ 
Lemma 5 For any M > 4, let a[^^^ be the coefficients defined by Eg. ( |T6l) . Then 
1. if M is even, then a^p ^'^ 7^ '^if-i-p /^'^ '^''^V P' 



2. if M is odd, then ap^'^^ 7^ '^M-p for o-ny p. 
roc 

ci'^M-2-p- expanding this polynomial we obtain that 



Proof. Since the polynomial defining the coefficients a^-^^'' is symmetric, we have ap^"^ 



AI-2 M-2-j , j 

j=0 k=0 ^ ^ i=0 

,=0 fc=0 t=o ^ ^ 



and so 



,{M) 



\M-2-k- 



p M~2-p 

E E ( 

fc=0 j=0 
p M~2~p 

E E ( 

fc=0 j=0 
p-l M-2-p 

E E (-1)"-^-"="' , , 

fc=0 j=0 

«rr'^ + af-'^ + (-i)''-' 



A; 



A;- 1 



k + j-l 
k 



sM-2 



s p M-3~p 



k=0 j=0 



+ (-1 



>M-2 
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Then the values of a^p^^ can be easily computed, see Table [H 



Table 1: Values of 



M 


"o 




4"> 


4"" 




4"> 




4 


1 


1 


1 











5 





1 


1 











6 


1 


2 


3 


2 


1 







7 





2 


4 


4 


2 







8 


1 


3 


7 


9 


7 


3 


1 



One may further observe that, for any even number M, a^^^'' 7^ cl^m-i-p ^p"^^ 

Therefore, the two parts in our lemma are indeed equivalent, hence we only need to focus 
on the case that M is even. 

For even M = 2K > 4 and 1 < p < M — 3, we always have Cg*^'' = 1 and 

AM) _ (M-1) (A/-1) 



a^'^^ = a';X'> + a\:''~'^ + l (27) 

ip_2 ^ + — 1 + Ctp. 



a%~'^ + aS"'^ - 1 + ai'!['^ + af -1 + 1 (2^ 



= 4^^-'^ + 2a(i^r'^ + af - 1. (29) 

We then show that {a^Q^\a^i^\--- ,a^^]} is an increasing sequence. It is obviously 
true for i^' = 3,4. Assume our claim is true for K = Kq, let's look into the case K = 
Ko + 1. We have 4'"^°+'^ = a^^^^'^ + 2af_'^"^ + 4'"^°^ - 1 for any 1 < p < Kq. Then 

r {2Ko+2) {2Ko+2) (2_ft'o+2)-| . , . . n- (2Xo+2) (2_ft'o+2) 

{Oq , a-[ ,■ ■■ ,a)^^_i \ IS also an mcreasmg sequence, bmce a)^^ — ciKo-i ~ 

(A2Ko)^ {2Ko) {2Ko)_.^_f(2Ko)^ {2Ko) {2Ko)_y. _ {2Ko) , (2Ko) _ {2Ko) _ {2Ko) _ 
\'^Ko-2^ '^"'Ko-l^^Ko ^1 l"/^o-3"'"'^"A'o-2"r"_fs:o-l ~ "-ft'o "'""/<o-l "/<o-2 "Xq-S ~ 

c^K^-i — (^^'k^~\ > 0' whole sequence {aQ^^"'*'^'', af^°^'^\ ■ ■ ■ , 0'k!^°~^'^^} is also an increasing 
sequence. This completes our proof. □ 



3.2 Universal subspaces of minimal dimensions 

It follows from Proposition [10] that any subspace of A^V which is spanned by less than 
M{M — 1)(M — 5)/6 basic trivectors eijk is not universal. It is natural to raise the question 
whether there exist universal subspaces spanned by exactly M{M — 1)(M — 5)/6 BSOVs. 
In the next proposition we show that this is indeed true when M is even, and verify it for 
odd M = 7,9,11. 

Proposition 6 Let V be a complex Hilhert space of dimension M and let e^jk = \i) A |j) A | A;) 
Jorl<i<i<k<M. 
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(a) If M = 2K > 6 is even, then the complex subspace spanned by all BSOVs Cijk except 
the K of them with indexes 

(1,3, 6), (1,4, 6) and (1, 2z - 3, 2z - 1) ?><i<K, (30) 

is a universal subspace of minimum dimension, M{M — 1)(M — 5)/6. 

(b) If M = 7, 9, 11, then the complex subspace spanned by all BSOVs Cijk except the M — 1 
of them with indexes 

(1,4,M),(2,5,M) and (z,z + 2,M) 1 < z < M - 3; (31) 

is a universal subspace of minimum dimension, M{M — 1)(M — 5)/6. 

Proof, (a) We can write an arbitrary \ip) G A^V as \ip) = J2i j k^ijk^ijk where the sum- 
mation is over all triples {i,j,k), l<i<j<k< M. First, by Theorem H] we may assume 
that the coefficients Cijk = whenever e^k is not a BSOV. For i < K, denote by Uj the 
subgroup of U(M) which fixes all basis vectors \j) with j ^ 2i — l,2i. All subsequent LU 
transformations will be performed by using these subgroups, and so the above mentioned 
property of the coefficients Cijk will be preserved. 

Second, we shall prove by induction that, by using only the Uj operations, we can achieve 
our goal to make Cijk = also when {i,j, k) is one of the triples (1, 3, 6), (1, 4, 6) or (1, 2i — 
3, 2i - 1) with 3<i<K. 

If K = 3, then M = 6. Note that the Hilbert space of 3-qubits can be isometrically 
embedded in this fermionic system as the subspace spanned by the BSOVs: 

\i) ® \j) ® |A;) -^\i + l) A |j + 3) A |A; + 5), i, j, k G {0, 1}. 

The desired result then follows from the well known fact [5] that by performing LU transfor- 
mations on any given pure state of three qubits one can make vanish any three coefficients 
(in the standard o.n. basis). This is because after the embedding, the three-qubit LU can 
be viewed as a special case of the fermionic LU, where we allow only the fermionic LU 
transformations given by block- diagonal unitary in U(6) with three 2x2 blocks in U(2). 

Now let A' > 3. By the induction hypothesis, we may assume that the coefficients Cyk = 
when {i,j, k) is one of the triples (1, 3, 6), (1, 4, 6) or (1, 2i — 3,2i — 1) with 3 < i < K. We 
can choose X G Uk such that X(ci,Af-3,A/-i|^ ~ 1) + ci^m-3,m\M)) oc |M). This means 
that the coefficient of ei,M-3,M-i in X ■ \ip) is 0. This completes the inductive proof. 

(b) This proof is computational. All computations were performed by using Singular 
PB] . We just have to apply ^ Theorem 4.2] as explained in the proof of Theorem HI The 
cases M = 7 and M = 9 were straightforward, and we omit the details. The values {q\p) 
that we computed for M = 7, 9, 11 are 48, 10368 and 12431232, respectively. In the case 
M = 11 we had first to eliminate the variable Xu in order to avoid the problem of running 
out of memory. We shall describe this elimination procedure in general. 

Thus, let M > 7 be an odd integer, and set M = 2K + 1. We define the polynomial 
ring P, its ideal X, and the inner product (■, ■) of polynomials as in the proof of Theorem 
m The polynomials go and p will be the same as in that proof, see Eqs. (fTTj) and (fT2|) . The 
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congruence of polynomials will be again modulo the ideal X. We shall also use the coefficients 
a^j^ , /c = 0, 1, . . . , M — 2, defined by the polynomial identity Eq. (ITB]) . (These coefficients 
are the same as the ones used in Lemma O) Other symbols that we are going to introduce, 
like /X, q and Si^k will have different meaning from the same symbols used in the proof of 
Theorem HI Since Cq^^^ = — 0; from Eqs. ( fTSj) and ( fT6l) we obtain that 

K M-3 
i=l k=l 

K K~l K 

i=l k=l f i=l 

where / runs through all functions {1, . . . , K} — > {1,...,/^' — 1} and 

_ fro- iTo l'^ ('fl^^^V^^~2-2/c , (M) M-2-2k\ 
Oi,k — l,"''2i-l"''2«j I U-fe ■^2i-l ^ ^M~2-k-^2i I 

= X2i-iX2is'ik, k = l,...,K-l. (33) 

Thus go = J2f1f where g/ = XiX2 - ■ ■ xm-iq'j and q'j = Y[f=i^if(i)- The characteristic 
polynomial of the subspace defined in part (b) is g = /igo where 

M-3 

H = (^Xi+X4 + xm){x2 + X5 + xm) JJ {xi + Xj+a + xm), (34) 

i=l 

see (l3Tll . Hence, this subspace will be universal if we can prove that (g|p) ^ 0. 

By using the expansion (fT2l) . we deduce that {q\p) = (yugo|p) = (Ai'gob) where n' is 
obtained from /i by setting xm = 0. After cancelling X1X2 ■ ■ ■ xm-i, we obtain {q\p) = {q'\p'), 
where g' = fi'q^, q'^ = J2fQ'f polynomial p' is defined by Eq. f|T2l) with M replaced 

by M — 1. Hence xm has been eliminated, it occurs in neither p' nor g'. For M = 11, we 
were able to compute {q'\p'). □ 

We remark that in the case M = 6 any subspace spanned by 5 BSOVs is universal. How- 
ever, in the case M = 7 the 6 BSOVs that we discard cannot be chosen arbitrarily, but there 
are several other good choices such as (1, 5, 7), (2, 5, 7), (1, 3, 7), (2, 4, 7), (3, 5, 7), (4, 6, 7). 

It should also be mentioned that, although it is not directly applicable, our work may shed 
light on the understanding of the N-representability problem [19], where the single particle 
eigenvalues are invariants under LU. For the = 3, M = 6, let us consider a concrete 
example of the universal subspace with dimension 5, spanned by 

6235; ^145, 6136, 6245, 6135, (35) 

where Cijk = \i) A |j) A\k). 

In other words, any three-fermion pure-state with six single particle states is unitary 
equivalent to 

\ip) = ae235 + beu5 + ccise + 0^6246 + ^6135, (36) 
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and the coefficients can be chosen as a,b,c,d > 0, z & C, and HV'IP = a^+&^+c^+c?^+|-2p = 1. 
Without loss of generahty, we can further assume that a > b > c. 

The one particle reduced density matrix pi of the state p := is given by (we choose 

the normalization Tr pi = 3) 



Pi = 



( b' + c' + 


2; P az 











\ 


az 


a^ + 




















+ + 


\z\^ bz 














bz* 


b^ + d^ 




















a^ + b'^ + \z\^ 


cz 


^ 











cz* 


c^ + d^ J 



We have pi = Ra (B Rb (B Rc, a direct sum of three 2x2 diagonal blocks. For x = a,b,c let 
Dx = det Rx- 

Since Tr Rx = 1 and Rx > 0, we have Dx G [0, 1/4] and the eigenvalues of Rx can be 
written as A^^ and 1 — A^^ with A^; = (l + Vl — 4Z}^.)/2 G [1/2, 1]. Let us denote the eigenvalues 
of pi arranged in decreasing order as Ai > A2 > ■ ■ ■ > A5 > Ag. Then Aj and A7_j are the 
eigenvalues of the same block Rx of pi. Thus the following result is a direct corollary of 
Proposition [6l 

Corollary 7 // the AjS are arranged in decreasing order, then Aj + Xj^i = 1 for i = 1, 2, 3. 

This hence gives an alternative proof for the N-representability equalities and inequalities 
in the case iV = 3,M = 6[20l|2ll[22l|23]. 

4 Universal Subspaces for > 3 

In this section we consider the generalization to > 3. Note that for even the total 
number of configurations is (^^) , and that the number of basic single occupancy states is 
2^(^). Consider the ratio r(iV,i^) = {^f^)/2^{^). We know that for any finite N,K, we 
have r > 1. If we fix A^, we will have limi^_^oo ''"(^7 -^'^) = 1- This means that the single 
occupancy states have measure 1 for large K, so one may hope that a similar result as for 
= 3 might hold for A^ > 3 at least when K is large. However we show that this is not the 
case. 

4.1 The SOV-subspace is not universal for > 3 

Proposition 8 Let M and N be integers such that M > 2N > 8. Let V be a complex 
Hilbert space of dimension M . Then there exist \ip) G /\^V which are not equivalent to any 

sov. 

Proof. Write M = 2K li M is even and M = 2i^ + 1 if M is odd. We fix an o.n. basis 
{\i) : I <i < M} oiV. Let Vi = span{|2? - 1), \2i)}, i = 1, . . . , i^, and let V^+i be the 
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1-dimensional subspace spanned by \M) if M is odd. Then the SOV-subspace, iS, of /\^V 
is given by 

S= Yl V,,AV,,A---AV,^. (37) 

l<ii<-<i]V<M 

Let G := GL{V) and let / : G x 5 ^ A^V be the map sending {A, ^ A ■ 
Let H be the subgroup of G which leaves invariant each of the subspaces Vi. Since S is 
if-invariant, we can form the algebraic homogeneous vector bundle G *hW with projection 
map p : G *h S — >■ G/H (see [271 section 4.8]). As the dimension of G *h S is D : = 
DimG + DimiS — Dim if, we have 

D = AK{K -l) + 2^(^y (Meven); 

D = 4i^' + 2^(^)+2^"^(^^l)^ (^odd). (38) 

Since / factorizes through the canonical map G x S G*h<S, we have Dim G ■ S < D. One 
can verify that D < ('J^), see Lemma O Thus G ■ S must be a proper subset of A^V. □ 

Lemma 9 For the dimension D defined by Eq. (155]) . we have D < 
Proof. We first prove the assertion for even M, i.e., that 

/(M, N) := - 2^ > AK{K - 1) (39) 

when M = 2K > 2N > 8. We use induction on A^. One can easily verify /(M, 4) > 0. 
Suppose /(M, A^) > holds for some A^ and that M > 2N + 2. The inductive step follows 
from the identity 

(A^ + 1) (/(M, A^ + 1) - /(M, A^)) = (M - 2A^ - 1)/(M, A^) + A^ ■ 2^ (^^^ . (40) 

which is easy to verify. Indeed, it implies that /(M, A^ + 1) > /(M, A^). 

Next we show the assertion for odd M = 2K + 1 > 2A^ and A^ > 4. For A^ = 4, we have 

,M\ . 4 fK\ , fK\ (2K + 1 

D-\ ] = AK^ + 2^ 



<^K, (^f ) + 2' (f j - Y ^ 2A-(3 - ) < 0. (41) 

where the first inequality is from Eq. fl39l) . For A^ > 4, we have 



Nj \NJ \N-1J \ N 

<-^-(C;)--(--))-((;",)-«(--))-r;' 

= AK{2 -K) <0, (42) 
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where the first inequality is from Eq. (139|) . So the assertion is true for odd M. This completes 
the proof. □ 
To conclude this subsection, we give the lower bound for the dimension of universal 
subspaces. 

Proposition 10 Let V be a complex Hilbert space of dimension M and : z = 1, 2, . . . , M} 

an o.n. basis of V . Let W be a complex subspace of A'^V spanned by some basis vectors 
\ii A 12 A ■ ■ ■ A In) , ii < 12 < ■ ■ ■ < iN- //Dim W < (^) — (^2^) then W is not universal (for 
the diagonal action of\J{M)). 

Proof. Let U(M) be the unitary group of V, identified with the group of unitary matrices 
of order M by using the basis {|^)}. Let / : U(M) x W — be the restriction of the 
action map U(M) x A^V — )■ A'^V. Since W is spanned by the basis vectors, it is invariant 
under the action of the maximal torus T of U(M) consisting of the diagonal unitary matrices. 
We can form the equivariant vector bundle U(M) *t VV. The map / factorizes through the 
canonical map U(M) x W — U(M) *t VV, and so we have 

Dim(U(M) ■ W) < Dim(U(M) *t W) = M(M - 1) + 2d, (43) 

where d is the complex dimension of W. As c? < (^) — (Y) , we obtain that Dim(U(M) ■ W) < 
2(5J) = DimR(A^y). Hence, U(M) ■ W must be a proper subset of A^V^, i.e., W is not 
universal. □ 
We do not know that whether the bound (^) — (*^) is sharp for universal spaces of 
general N > 2, M. For = 3 and M even, the bound is sharp by Proposition [6] (a). 



4.2 BCS states are not LU-equivalent to single occupancy states 

We construct a concrete example that is not LU-equivalent to any SOV, i.e. the BCS 
states. The BCS states are named after Bardeen, Cooper and Schrieffer for there 1957 
theory using the states to build microscopic theory of superconductivity [28]. Here we 
consider the finite dimensional version of the BCS state \iPn,m), which is an iV-fermion 
state, i.e. \4'n,m) ^ A^V with Dimy = M, where both and M are even (see e.g. [29]). 
The M states are 'paired' in a sense: say as {1, 2}, {3, 4}, . . . , {M — 1, M}. In each pair, the 
two states are either both 'occupied' by a pair of fermions or both 'empty'. The BCS state is 
an equal weight superposition of all those 'paired' states. Intuitively, BCS states are 'always 
double occupancy', so that they should be among the 'hardest' ones to be transformed by 
LU to some single occupancy state. 
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We list a few BCS states: 



K 

|V'2,m) = J]|2z-1) A|2z), 

|V'4,m) = - 1) ^ ^ - 1) A |2z2), (44) 

|^7V,m) = |2zi-l)A|22i)A---A|22^/2-l)A|22jv/2), (45) 

«l<«2<---<«JV/2 

where the G {1, . . . , iC}. 

We shall use exterior multiplication and partial inner products with the BCS states, e.g., 
\a) A \iIjn,m) {0'\'^n,m) where \a) G V. To simplify this computation, let Ui,...,Uk be 
2x2 unitary matrices with det Ui = 1. So U = (B^iUi is a M x M unitary matrix. One 
can verify that 

A^ t/|2i - 1) A |2i) = \2i-l)A\2i), (46) 
where i = 1, . . . , K. Hence we obtain the stabilizer formulas 

A^'UlijNM) = \i^N,M)- (47) 

Consequently, in the expressions \a) A \iPn.m) and {a\ipj\f^M), after an LU transformation we 
may assume that \a) is a linear combination of the \2i — 1), i = 1, . . . ,K. 

By Lemma ^ ii N = 2 then any BCS state is LU-equivalent to some SOV. In Theorem 
[T2] we will show that this is not the case for > 4. For this purpose we give a general 
criterion of deciding whether a state is LU-equivalent to an SOV. We shall denote by pu the 
bipartite reduced density matrix of p. 

Lemma 11 Let {ip) G A'^V, Diml^ = M > 2N , he an antisymmetric state and let 
p = Then is LU-equivalent to an SOV if and only if there is an o. n. basis 

|ai), . . . , |aAf) of V such that pi2(|a2i-i) A \a2i)) = for all i with 2i < M . 

Proof. Sufficiency. By replacing with an LU-equivalent state, we may assume that 
\ai) = \i) for all i. We can write 

1^) = A--- Aijv). (48) 

ji<---<ijv 

If ii = 2s — 1 and 12 = 2s, the hypothesis implies that |2s — 1) A |2s) G ker Tr3^...^Ar = 
kerpi2. So the coefficient ... = 0. Similarly, Ci^^...^j^ = if {2s — 1,2s} C {i^, . . . ,ij^} 
for some s. So is a single occupancy state. 

Necessity. Suppose \ip) is LU-equivalent to a single occupancy state \(p), i.e., \ip) = 
A^U\(p) with a unitary U. By definition of SOV we have ((/?|(|2s — 1) A |2s)) = 0. This is 
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equivalent to ((2s — 1| A {2s\)a{\2s — 1) A |2s)) = 0, where cr = \(p){(p\. By tracing out all 
but the first two systems, we obtain ((2s — 1| A (2s|)(Ti2(|2s — 1) A |2s)) = 0. It follows that 
|2s — 1) A |2s) G ker(Ti2. Thus, the assertion holds with |aj) = U\i) for all i. This completes 
the proof. □ 
Note that the proof can be easily extended to the case in which LU is replaced by the 
diagonal action a'^V where V is invertible. Now we prove the main result on BCS states. 

Theorem 12 Let M = 2K and N be even integers with K > N > 4. Then the BCS state 
\'4^n,m) is not LU-equivalent to any SOV. 

Proof. We start with the case iV = 4. Suppose |^/'4,m) is LU-equivalent to an SOV. By 
Lemma [TH there is a nonzero decomposable 2- vector \ip) = \a) /\ \b) such that the partial 
inner product (v?|V'4,m) = 0. By the simplification mentioned beneath Eq. ( 14 7p . we may 
assume that \a) = XljLi '^2i-i|2j — 1), ai = 1, and \b) = X]jl2^ib)- ^ special case of 
formula ([2]) we have 

{aAb\viAv2Av3AVi)= ^ {-iy+^-\a A b\vi A Vj)\vk A Vi) , (49) 

l<j<j<4 

where {i,j,k,l} = {1,2,3,4} and k < I. Recall from subsection 12.31 that i = {i,i — l} when 
i is even and i = {i,i + Ijwhen i is odd. For convenience, we shall write i < j ii i < j and 
i ^ j. By using the formulae in Eqs. ( H9l) and (Hll), we conclude that 

{aAb\iA]) = Q if 2 < j, (50) 

where j G {1, . . . , M}. By setting z = 1 in Eq. f l50|l . we conclude that bj = for j > 2 
and so \b) oc |2). Next, by setting z = 2 in Eq. f l5U]) . we conclude that also aj = for j > 1 
and so \a Ab) oc |1 A 2). As (1 A 2|?/'4_Af) 7^ 0, we have a contradiction. Thus |V^4,m) is not 
LU-equivalent to an SOV. 

Next, we use the induction on to show that there is no nonzero decomposable 2- vector 
\ip) = \a Ab) such that the partial inner product {(p\4'n,m) = 0. By Lemma [TTl this implies 
that the BCS state \iI>n,m) is not LU-equivalent to any SOV. 

Since we have already verified the first case, = 4, let us assume that the assertion holds 
for \'iPn-2,m) with A^ > 6 and M > 2{N — 2). Suppose there is a nonzero decomposable 
2-vector \lp) = \a Ab) such that {^p\ipN,Ai) = with M > 2N. Let 

M M 

I") = ■=\(^) + ^J'l-^'^' ^^^^ 

j=l j=M-l 
M M 

\h) = Y.^3\j) ■■= \P) + E ^^-i^')- (52) 

3=1 j=M-l 

Since {(p\iPn,m) = 0, we have 



= (M-1|A(M|(((p|^^,,,)) = (^|(^((M-1|A(M|)|V^^,m) 

= (v5|^^Ar-2,M-2) = (« A /3|^Ar_2,A/-2)- (53) 
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As M — 2 > 2{N — 2), the induction hypothesis imphes that |a A /3) = 0. Thus \a) and 
1/3) are hnearly dependent. Consequently, we may assume that = 0, and moreover that 
\b) = \M) and that all = 0. Then by expanding the left hand side of (a A b\ip]\f^M) = 0, 
we obtain a contradiction similarly as in case = 4. So there is no nonzero decomposable 
2-vector \ip) such that {(p\'iI>n,m) = 0. This completes the proof by induction. □ 

5 Summary and Discussion 

We have discussed universal subspaces for local unitary groups of fermionic systems. We 
have shown that for = 3, the SOV-subspace is universal. Furthermore, for M even, we can 
always find a universal subspace, contained in the SOV-subspace, whose dimension is equal 
to the lower bound M{M — 1)(M — 5)/6. Although our main tool is a natural application of 
Theorem 4.2 in [6], which can be used in small dimensions to construct universal subspaces 
by computers, the analytical proof we obtained for the general case is far from trivial. In 
fact, some special features of polynomials are used, which do not generalize to settle the odd 
M case. 

We have also shown that, for A^ > 3, not all fermionic states are LU-equivalent to a single 
occupancy state. Our argument is based on dimension counting. For M even, we give BCS 
states as concrete examples that are not LU equivalent to any single occupancy state. This 
is intuitive as BCS states are always paired so they are most unlikely to be transformed with 
LU to something unpaired. Given that for fixed A^ almost all states are single occupancy in 
the large M limit, the BCS states are also among the 'measure zero' states. 

We wish our results shed light on further study of entanglement properties on fermionic 
system, as well as other properties of fermionic states such as the N-representability problem 
(as we discussed for the A^ = 3, M = 6 case). We also leave some open questions. One of 
them is to understand the achievability of the dimension lower bound for A^ = 3 and M odd, 
or for even the case of general N > 2, M, which is worth further investigation. 
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